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Abstract 

Quantum fluctuations in a current biased Josephson junction, described in 

terms of the RCSJ-model, are considered. The fluctuations of the voltage 

and phase across the junction are assumed to be initiated by equilibrium 

current fluctuations in the shunting resistor. This corresponds to low enough 

temperatures, when fluctuations of the normal current in the junction itself 

can be neglected. We used the quantum Langevin equation in terms of random 

variables related to the limit cycle of the nonlinear Josephson oscillator. This 

allows to go beyond the perturbation theory and calculate the widths of the 

Josephson radiation lines. 
PACS numbers: 74.40. +k 
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INTRODUCTION 



Fluctuations in an equilibrium system are considered to be "classical" if the fluctuation 
frequency is lower than the the system temperature, %u <C T. For higher frequencies the 
fluctuations are "quantum" and contain a contribution called zero-point fluctuations (ZPFs). 
Due to the ZPFs there are fluctuations even at zero temperature, when the system is in its 
ground state [|TJ. 

The intriguing question: can these zero-point fluctuations be detected, has recently at- 
tracted much attention. The answer depends on what is the detector and what means 
detection. Measuring Casimir forces between two bodies, both at T = 0, is definitely a way 
to detect ZPFs of the electromagnetic field surrounding the bodies. On the other hand there 
is no radiation from a body at T = and an antenna placed nearby will receive no signal, 
which means that current ZPFs in the body do not radiate. 

The situation is even more complicated if the system is not in equilibrium. Probably the 
simplest example is a gas of excited atoms in a cavity with no radiation in it (the temperature 
of the radiation is zero, while the temperature of the gas is not). When de-excited, atoms 
emit spontaneously cavity photons, which is in fact interaction with ZPFs of the cavity field. 
The de-excitation of the atom can be considered as a signal of the presence of the ZPFs in 
the cavity. 

This paper is stimulated by experiments with nonequilibrium pumped Josephson junction 
systems P||}4|], where ZPFs play an important role at low temperatures. In ZPFs were 
detected in a dc-current biased Josephson junction, measuring the low frequency noise of 
the voltage V(t) across the junction. The voltage noise was generated by intrinsic current 
fluctuations SI(t) in a shunting resistor R. The frequency of the measured voltage noise (of 
the order 100 kHz) is much below the temperature T of the resistor (1.6 K and 4.2 K) but 
quantum phenomena are important because this low frequency is mixed due to nonlinear 
effects with the Josephson frequency uoj (between 10 and 500 GHz), which is comparable or 
even higher than the resistor temperature. 
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The theory of quantum fluctuations in a dc-current biased Josephson junction was given 
in ||, based on the Langevin equation for the phase difference across the junction ip(t), 
where the random forces are intrinsic current fluctuations 51(f) in the resistor. We decided 
to revisit the problem because of following reasons. In || the authors used the theory of 
classical fluctuations developed in simply replacing in the classical Langevin equation the 
classical spectral density of the random forces, (dl 2 )^ = T/tcR, by its quantum equivalent, 
(huj/ZirR) coth.(hu/2T). This quantum spectral density, containing ZPFs, is generated by 
the symmetrized correlator ((l/2)[8I(t')5I(t) + SI(t)8I(t')]), which is, however, not relevant 
in the problem considered. This is because using the symmetrized correlator for the random 
forces in the Langevin equation imply that all the correlators calculated from this equation 
are also symmetrized. On the other hand, the symmetrized correlator of the voltage V(t) 
do not represent the measured voltage noise. 

This problem was addressed in [[F],|]|§. It was shown that if the device, measuring current 
noise, is a "resonator" at zero temperature, with resonant frequency u, the measured signal 
is proportional to the Fourier component S(u) of the nonsymmetrized current correlator 
(j(O)j(t)), with uj > and the dc component in the current j(t) being subtracted. S(uj) 
is real, but not symmetric in uj, and hence the condition uj > is essential. If the current 
carrying system is in equilibrium at zero temperature, one finds S(u) = for uj > 0, which 
means there is no signal created by ZPFs. It is because the signal S(u) is in fact proportional 
to the power spectral density spontaneously radiated by the system. The situation is similar 
in quantum optics [|l(|, where ZPFs can create radiation only in nonequilibrium systems, 



like parametric amplifiers. 

The second reason to revisit the theory is as follows. In || and || the fluctuations were 
calculated using perturbation theory for the phase (p(t) , i.e. assuming that (Sip 2 )^ is small 
and proportional to (SI 2 )^. This theory diverges near uj = and the harmonics of Josephson 
frequency u = kuj, k = 1, 2, .., which are exactly the points of interest. 

To account for all the above mentioned circumstances we use the quantum Langevin 
equation for ip(t), formulating it according to the general recipe given in [fil l. With no 



random forces the solution of the Langevin equation for the phase, Tp{t) = f{ujjt + aco), 
contains an arbitrary initial phase a , since the equation for Tp{t) is invariant with respect 
to time shift. Due to random forces this initial phase acquires a fluctuating contribution 
5a(t), and the Langevin equation can be rewritten for 5a(t). Using this approach we avoid 
the perturbation theory for tp(t), which allows to calculate the shape of the radiation lines 
at the Josephson frequency harmonics kuoj. 

The width of the Josephson emission lines was considered in the well known papers Jl2| , 
however a white spectrum was assumed for the random forces. This assumption is employed 
in almost all papers dealing with noise influence on Shapiro steps |[L3|| , on the impedance 
of microwave driven junctions and other phenomena in Josephson junctions (more 

references are given in |fT6|JT7f| ) . There are only few exceptions. In [p~8] ? p~9|1 the white noise 
with correlator ~ 5(t) was replaced by a dichotomous telegraph noise with an exponential 
correlator ~ exp(— yt). In p(J the sum of these two was considered. Full consideration for 



the quantum correlations, considering the rounding of the I — V- curve, was given in [pll] . 



I. QUANTUM LANGEVIN EQUATION 

The equations describing a current biased Josephson junction in the RCSJ-model are as 
follows 

. • ,s V(t) „dV(t) T er/ x d<p(t) 2 (3 . « 

I c sin tp(t) + + = I P + 5I(t), -^2 = -V{t). (1) 

Here I c and I p are the critical and bias currents, R and C are the resistor and capacitor of 
the model, ip(t) and V(t) are the phase and voltage difference across the junction and SI(t) 
are the intrinsic current fluctuations. The phase and voltage differences contain fluctuations 
generated by SI(t) and in what follows we write <p(t) = Tp(t) + 5ip(t) and V(t) = V(t) + SV(t), 
separating explicitly the fluctuating contributions. 

We neglect the quasiparticle current in the junction, assuming low enough temperatures, 
which means in terms of the model that R is just the external shunting resistor. With this 
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assumption SI(t) are the equilibrium current fluctuations in the resistor R, the only intrinsic 
fluctuation source in the RCSJ-model This simplifies the problem enormously, since 

the current fluctuations in the junction are state dependent (they depend on Tp(t)) and hence 
are non-stationary for time-dependent voltage V(t) p3|j2^ , |25[] . 

We eliminate from Eqs.([L]) the potential V, multiply it by h/2e and introduce E c = 
(H/2e)I c . As a result we have 



h \ 2 d 2 (p l / h\ 2 dip % h 

— — — r H — — ; — h E c sin ip I v = — . 

2e dt 2 R \ 2e ) dt * 2e p 2e 



C \-) -^ + ^(-1 -JI + E c^<p--I p = -6I. (2) 



This equation is a standard form of a Langevin equation for the phase (p as the particle 
"coordinate" and a "potential" having the dimension of energy, the current fluctuations 
being the random forces. Comparing with the quantum Langevin equation in [|lTJ we find 
the correlator 

(8I(t')8I(t)) = 4; r dujnuJ N^e-^*-*') + (N(u) + l)^^'^] (3) 
ttR Jo 1 J 

and the commutator 

[SI{t'),SI{t)} = =-\ duhuu sin u(t-t') = j-5(t-t'), (4) 
ttR Jo R dt 

where N(ou) = [exp(hu /T) — . It is important that 8I(t) is a Gaussian process (in the 

quantum case one has to preserve the order of operators) and the commutator Eq.([|) is 

a c-number. (When the thermal bath responsible for the random forces is chosen to be a 

collection of harmonic oscillators JT1J the Gaussian properties of the random forces follow 

from the Gaussian properties of a harmonic oscillator in equilibrium. However, probably, 

this is a more general property, since many degrees of freedom of the thermal bath contribute 

to the random force.) 

We will use the following notations for Fourier components and spectra of correlators 

/+oo 
due-^A m (A% = (AJt = Al ( 5 ) 
-oo 

f-\-oo 

(A(t')iA(t)) = / dcoe-^-t'^AX, (4,^) = S(u - ^)(A\, (6) 



where t means "hermitian conjugate" . In these notations the spectrum of equilibrium current 
fluctuations is 

= ^N(u). (7) 

The width of the spectrum is T/% and the corresponding correlation time is fi/T. At T = 
this spectrum has no Fourier components with u > 0, which means that energy can not be 
extracted from the thermal bath. The spectral density of the symmetrized correlator, which 
contains ZPFs, 
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does not have this property. In the classical case, when T ^> Too and \t — t'\ ^> h/T, the 
random forces have a "white" spectrum, i.e. 

T 9T 
(51% = — , (SI(t')5I(t)) = —5(t - 1>). (9) 



II. GENERAL CONSIDERATIONS 

We recall first the properties of the junction when one neglects fluctuations. The solution 
of Eqs.(|l|) with 51 = differs qualitatively in cases when the bias current I p is smaller or 
greater than the critical current I c . We will consider only the case of the non-stationary 
Josephson effect, when I p > I c . In this case the voltage V(t) and the phase Tp{t) are periodic 
functions of time (the last one modulo 2n), the period defining the Josephson frequency ujj. 
These functions can be presented as follows 

Tp(t) = f(x), V(t) = I C R g(x), x = ujt + oq, (10) 

where f(x) and g(x) are non-dimensional functions 

g(x) = ^ 9k e- lkx , f{x) = £ f k e- ikx + x, (11) 

k k 
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ao is an arbitrary initial phase, appearing because Eqs. ([!]) are invariant with respect to time 
shift. Eqs.flTOD define the limit cycle of the nonlinear Josephson oscillator, i.e. the trajectory 
of the system in the phase space (<p,V). The time-average voltage across the junction is 
V = I c Rgo- 

As was explained in the Introduction the measured voltage noise (proportional to the 
radiated energy) is 

1 r+oo 

S{u>) = - j ^ dte^(v(0)v(t)) (12) 

with v(t) = V(t) — Vq and u > 0. When fluctuations are neglected, v(t) is a non-random 
periodic function. Considering the initial phase ao to be random within the interval (0, 2ir), 
we convert the periodic function v(t) into its random equivalent, i.e. a stationary random 
process with an equidistant discrete spectrum. Using 



exp [-i(k - k')a ] = 5 k ,k', (13) 
we average v(t')v(t) over ao and obtain the relevant correlator 

(v(t')v(t)) = (I C R) 2 E \gk\ 2 exp[-ikuj(t - t')] (14) 
and the noise spectrum 



v 
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• S(u) = {I c Rf E \9k\ 2 5{uj - kuj), lu>0. (15) 

We turn now to the effect of fluctuations, when the random forces SI(t) tend to move 
the system away from the limit cycle and destroy the periodicity of the system dynamics. 
We assume the random forces to be weak. If the system is shifted by the random force 
perpendicular to the cycle, a "restoring force" exist, which tends to bring the system back 
to the cycle. However no such force exist when the system is shifted along the cycle, since 
this corresponds to a change of arbitrary initial phase «o- To account for these two effects 
we look for a fluctuating solution of Eqs.(|l]) in the form 

<p(t) = f(wjt + aa + 6a(t)), V(t) = I c R[g{ujt + a + Sa(t)) + 5£(t)], (16) 



with two unknown random functions 5a(t) and 5£(t) (which are operators in the quantum 
case). The last one describes the shift of the system perpendicular to the limit cycle and since 
in this case a restoring force acts against the random forces, <5£ (t) is small and proportional to 
5I(t); it contains the same frequencies as 51 (t). Contrary, 5a(t) can grow unlimitedly even for 
small SI(t), if the random force contains the resonance frequency ujj or its harmonics. Hence, 
8a(t) = 5a(t) + 5a(t), where Sa(t) is proportional to SI(t) and contains the frequencies of the 
random force, while 5a(t) is secular and slow, i.e. contains only low frequencies proportional 
to the amplitude of the random force. The function Sa(t) describes the so called diffusion 
of the initial phase and its derivative describes Josephson frequency fluctuations, which are 
responsible for the broadening of the Josephson radiation lines, 

Sujj{t) = ^-5a{t) = 5a(t), 5a(t) = f dt'5a(t'). (17) 

In terms of the new variables the voltage noise is 

(v(t')v(t)) = (I C R) 2 x (18) 



£ g* v exp[ik'(ujjt> + « + 5a(t'))} + 5£(t') 

k'^0 



9k exp[-ik(ujt + a + Sa(t))] + 6£(t) 

kjtO 



As we will see later {Sa(t), S£(t)} is a Gaussian process and to find the voltage noise one 
has to calculate the correlators (da(t')da(t)), (dat(t')5£(t)) and (5£(t')5£,(t)). 

The effect of the random forces on the noise spectra Eq . (|12D is twofold: the monochro- 
matic Josephson lines in Eq. fll~5|) acquire a width, and a background for all frequencies uj > 
appears. Weak random forces means that the width of the Josephson lines is small compared 
to the distance between them and that the background contribute to the integrated noise 
less than the lines do. 

Of special interest is the background 5(0) at low frequencies lo ujj. To calculate the 
background noise at frequencies far enough from the resonant ones to = kutj, k = 1, 2, ... one 
can employ perturbation theory assuming Sa(t) to be small, while to find the width of these 
lines a more sophisticated approach has to be used. 
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III. LANGEVIN EQUATIONS FOR AN OVERDAMPED JUNCTION 



Substituting Eqs.fllO") into Eqs.fll]) with 51 = we find the relations between the functions 
/ and g as follows 

g(x) + sin f(x) + ujrg'(x) = p, uJjf(x) = u Q g(x), (19) 

where p = I p /I c , &o = (2e/ft)I c i£, and r = RC. To obtain the Langevin equations for 5£(t) 
and 8a(t) we substitute Eqs.(^) into Eqs.fll]) and use Eqs.fll9"|). As a result we find 

r5i(t) + 5£(t) + rg'(x)8d(t) = 8I(t)/I c , f'(x)8a(t) = w 8£(t) (20) 

with x = ujjt + ao + 5a(t). (Some caution is needed in the quantum case since 5a and 5a 
do not commute. However the main (secular) part of 8a is a slow function of t and hence 
8a(t) and 8a(t') effectively commute when t' = t.) 

In what follows we will consider the overdamped Josephson junction, assuming the ca- 
pacitance C to be small, i.e. ujqt <C 1. In this case we can skip in the first of Eqs. (|20|) the 
term with 8a(t), since, according to the second of these equations, it is a small correction to 
the term 8£(t). We can also use the functions f(x) and g(x) as for a junction with C — 0, 
in which case p3[ the Josephson frequency is uj = co a with a = (p 2 — l) 1 / 2 and 



g(x) = (p 2 -l)/(p + sin x); g k = {-i) k a{p-a) k ) k > 0; g„ k = g* k . (21) 

The functions f(x) and g(x) for an overdamped junction are related as follows 

<rf'{x) = g(x), g{x) + sin f(x) - p = 0. (22) 

As a result we have the following Langevin equations 

r8i(t) + 6£{t) = 8a(t) = —[p + sm(ujt + a + 5a{t))} 6£{t). (23) 

l c a 

From the first of these equations we find immediately 

i c 1 — 2C<JT i c 1 + (CJTJ^ 



where 5 J are current fluctuations shunted by the capacitor. As mentioned already 8a(t) = 
5a(t) + 5a(t), where 5a (t) is fast oscillating but small, while 5a (t) is slow but generally 
not small. When performing the Fourier transform of 5a(t) in the second of Eqs.(p3|) one 
can neglect in the right hand side the small 5a(t) and replace the slow 5a(t) by a constant, 
absorbed in a®. With this approximation 



5a tl 



-j- I -5 Jui + 7^-fexp[ia ]5J w+a;j - exp[-iaro]<yj w - W jl } • (25) 
l c la lia L J J 



If one takes into account the slow variation of 5a(t), then 5J W is replaced by its convolution 
with the narrow spectra of 5a(t), but this effect can be neglected for weak enough random 
forces. In the time domain Eq. (|25|) corresponds to 

5a(t) = — [p + sin(uj jt + a )}^^-. (26) 
a l c 

To calculate {5a 2 ) UJ we average over ao using Eq. (|13D and obtain 

(5a\ =(y) 2 { (a) 2 {5J2) " + 4^ + ( Sj2 ^->] } ■ ( 27 ) 

It is convenient to present 

■77pr= E he~ ikx , b = P -, b ±1 = ±— ; b k = 0, (|*| > 1). (28) 

J \^ ') k=—oo ® ^ 

Using this representation and the intrinsic current fluctuation spectrum given by Eq.(^) we 
have 

where g = {h/e 2 )R~ l is the non-dimensional conductance of the junction resistor. It what 
follows we assume g > 1, and we will see that this is the main condition for the fluctuations 
to be weak. We will not consider very strong pumping and assume p ~ 1. Then it follows 
from ujqt <C 1 that also ujjt <C 1 and one can simplify Eq. (p9|) , 



V* 2 )<* = 4 n( ^ } S N 2 {^ML u - kuj > n( w ) = 1 + ( ^ r)2 , (30) 

where n(u;) is the shunting factor. It is obvious from Eq.fl2"T|) that approaching the threshold 
p = 1 of the non-stationary Josephson effect, when a — > 0, the fluctuations increase. 
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IV. PERTURBATION THEORY 



We consider shortly the situation when Sa(t) can be assumed to be small. As we will see 
no problems appear in perturbation theory if one neglects the capacitance of the junction 
and put r = 0. Substituting in the second of Eqs. (|T6|) the Fourier expansion of g from 
Eqs.(|TT|) and expending in Sa(t) we find the voltage fluctuations 

5V(t) = V(t)-V(t) = I c RY,{-ik)g k e^[-ik(ujt + a Q )]5a{t) + R5I{t). (31) 

k 

Using 5a^ = 5a w /(—iuj) we have for its Fourier transform 

SV U = I C R £ kg k eM-ika } + RSI U . (32) 

k UJ — KLOj 



We substitute here Eq . (|25|) and shift the summation over k, obtaining 

SV U = ^exp[-ika ]Z k (u)6I u -. kuj , (33) 
k 



where the impedance of the junction is [|] 

v ( \ o I x , P kg k 1 

Z k (u) = R < Ofe o H : h 77— 

I a uj — kujj lio 



(k + l)g k +i {k - l)g k -i 



(34) 



uj — [k + l)u j 00 — (k — 1)lo j 

Averaging over the initial phase ao we have the general result of the perturbation theory as 
follows 

(WX = E \ZkW\ 2 ($I%-ku,j, lo>0. (35) 
k 

In the classical case one can see, using Eq.((|), that the above result agrees with ||. However 
in the quantum case the voltage noise is not given by the symmetrized current spectral 
density, as it was obtained in 0. It follows from Eq . (|32|) that perturbation theory diverges 
at u) — kujj for k — 1,2, i.e. at the frequencies of the Josephson lines. This is because, 
as already mentioned, at these frequencies 5a(t) is not small. Note, that in 0, where 
perturbation theory was developed for 5(p(t), there was also an artifact divergency for (Sip 2 )^ 
at uj = 0. 
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Using g k from Eq.(|2~l"D one can check that Zk(0) = when \k\ > 1 and find the low 
frequency uj <^ujj voltage noise to be 



S(0) = (5V 2 ) = \Z (0)\ 2 (8I 2 ) + ^(O)! 2 ^ 



(36) 



with \Z (0)\ 2 = R 2 (p/a) 2 } |^i(0)| 2 = R 2 /4a 2 , giving 



5(0) 



RTp 2 



7r a c 



(37) 



2p 2 2T V 2T 

In this special case the symmetrized current spectral noise enters, and this is why our result 
for 5(0) agrees with that obtained in B. 



V. DEPHASING FACTOR OF JOSEPHSON LINES 

To find the width of the Josephson lines one has to go beyond the perturbation theory 
to eliminate the singularities of (SV 2 )^ near the Josephson frequency and its harmonics. If 
5Vuj is calculated from Eq.(j3^) in the vicinity of kuj, (k = 1,2...), the term RSI^ can be 
neglected, and (SV 2 )^ is a double sum over k and k' . The strongest singularities (u — kujj)~ 2 
emerge from the diagonal terms with k! = k, while the non-diagonal terms create weaker 
singularities (u — kujj)' 1 . Comparing the diagonal and non-diagonal terms reveal that when 
\uj — ku)j\ <C ujj the diagonal terms dominate. Hence, to find the width of the Josephson 
lines kujj on can neglect in Eq.fllSD the terms 5£(t) and 5£(t') and pick-up in the double sum 
over k and k! the terms with k! = k only, considering 

(v(t')v(t)) k = (I c R) 2 \g k \ 2 exp[-ikujj(t - t')] (exp[ik5a(t')}exp[-ik5a(t)]) . (38) 

As a result the delta-function in the spectrum Eq.(|l|) is replaced by a form-factor 

poo (It 

M")= —exp[iut)D k (t) + c.c, (39) 
Jo 2ir 

where v = uj — kujj is the detuning from the center of the line and the dephasing factor is 

D k (t - t') = (ex-p[ikSa(t')) exp[-ik5a(t)]) = D k {t' - t)*. (40) 
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The symmetry relation for the form-factor follows from the invariance of the average with 
respect to time shift. Since, as assumed, the fluctuations SI(t) are small, the dephasing 
factors decay slowly with t. 

As follows from Eq.(^6|) and the properties of 81 \t) described in sec.|, that 8a(t) and 
5a(t) are Gaussian processes and their commutators are c-numbers. Using the well known 
identity e A a B = e A + B e \i A > B \ one can write the operator product entering the dephasing 



factor Eq. fl40|) as 

exp[ik8a(0)} exp[— ik8a(t)} = exp[ik(8a(0) — 8a(t))} exp 



l -k 2 [8a(0),8a{t)} 



(41) 



Using the Gaussian properties we calculate the averages and have for the dephasing factor 



D k (t)=exp k 2 C(t) exp -k 2 K(t) 



(42) 



where 



C{t) = hsa(0),5a(t)] = ]- f dt x f dt 2 [5a{t l ),8a{t 2 )}, 



(43) 



K(t)= 1 -([5a(0)-5a(t)] 2 ) = ± 



dt'8a{t') 



n 2\ 



l f dtx f dt 2 {5a{t 1 )5a{t 2 )). 
I Jo Jo 



Presenting the commutator entering C(t) as its average (8a(ti)8a(t2) — 8a(t 2 )8a(ti)) , and 
substituting the Fourier representation of the correlator (6a(ti)6a(t2)) , we find 



K(t) = dco g(u) 



1 — cos cut 



CO 



C(t) = i du q(u) 



sin uot 



CO 



with 



(44) 



(8a 2 ) u + {56? 



q(uj) = - - (8a 2 )- 



(45) 



Using Eq . (|30|) and the relation (— cj)N( — uj) = u(N(u) + 1) one can check that 

4 



i^ 2 U = — n(«) £ \h\ 2 MN(u>) + i]} w ^ u _ ^ 

n 9 k 



and result, 



(46) 



g(u) = — H(«) £ \b k \ 2 {u [2N(u) + 1]} U ^ u _ kuj , 

n y k 



(47) 
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A A On 2 4- 1 

q(u) = -U{u) £ \h\ 2 M w => u -ku,j = -^^- n H- («) 

q{oj) contains only ZPFs and is temperature independent. This is not the case for g(u>), 
which for high and low temperatures is, correspondingly, (for u > 0) 

4 2p 2 + 1 . . , / fiuj\ , , 

,M = -^nH.coth(^-j, (49) 
= — -Ln(w)[e(wj - -u) + (2 P 2 + 

where Q(ou) is the step function. Because of the ZPFs g(uS) ^ at T = 0. The shunting 
cut-off is not crucial in calculating C(t); this integral converge even if one replace II (u) => 1. 
Contrary, the integral K(t) converges only due to this cut-off. 

VI. QUANTUM DYNAMICAL NARROWING 

As can be seen from the previous section the form-factor of the Josephson line is given 
as an integral of the type 

$(!/)=/ — exp[iut]D(t) + c.c. = —coB[ut + H(t)]exp[-W(t)], (50) 
Jo 2ir Jon' 

where the dephasing factor is D(t) = exp[— W(t) + iH(t)]. The functions W(t) and H(t) are 
given by their spectral representations 

W(t) = I™ dujG{uj) l ~ C0 2 SU \ H(t) = f°^QM^, (51) 

JO LO z JO UJ Z 

where both spectral densities G(u), Q(u) contain a small pre-factor (due to g~ l ) and 
G(uj) > 0, while sign Q(u) = sign u>, Q(0) = 0. This is the quantum version of spectral 
line dynamical narrowing. (In the classical theory H{t) = 0, hence the line is symmetric 
in v, and G'(0) = 0. One can see from Eqs. (fj9|) that this property indeed holds for high 
temperatures but is not valid for zero temperature). 

The large time asymptotic of W(t) is obtained replacing G{uS) by G(0), giving 

W(t) = yt, 7 = |G(0). (52) 
14 



Since Q(0) = the function H{t) is finite at t — > oo. If one replace W(t) by its large time 
asymptotic and neglects H(t), the form factor becomes a Lorenzian 

m = --r—^- (53) 

As we will see, this line shape is valid only for small enough v. To get an expression valid 
also for large v and to estimate the non-Lorenzian contributions we proceed as follows. We 
present 

W(t)= 7 t + W(t), W(t)= du[G(w)-G(0)] = , (54) 

JO U) z 

and expand exp —W(t) = l — W(t). Neglecting second and higher order terms in W(t) and 
H(t), we present the form-factor as a sum of a symmetric and antisymmetric contributions 

$(!/)= $»+$ a (l/) (55) 



with 



$ » = ~ / -exp[- 7 t]cosz/t W(t), (56) 

7T V 1 + Y Jo 7T 

roo 

$ a (z/) = - / _ exp[-7f] sin vt Hit). (57) 

Jo 7T 



To calculate the non-Lorenzian parts of the form-factor introduce functions 

2 r°° 

A c (lj,is) = / dt cos vt exp [— 7 t] ( 1 — cos ut) = (58) 

n Jo 

2 7 ic> 2 (3z/ 2 - u 2 - 7 2 ) 

n ( 7 2 + z/ 2 )( 7 2 + (u - cu) 2 )( 7 2 + (z/ + cu) 2 ) ' 

2 Z" 00 2 l^yuiv 

A 8 (u,v) = - rftsinz/texp[- 7 t]sinwt = --— — — ■ — — . (59) 

TV Jo 7T (7^ + (i/ — o;) 2 )(7 2 + (y + cj)" 1 ) 

With this definitions we have 

= + I / o °° ^[GH - G(0)]A c (u,, „). (60) 

*.(") = -«/ —^Q(oj)A s (u,u). (61) 
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For small 7 in the vicinity of uo = v > both A CjS (uj, v) are smeared ^-functions, 

A CjS (o;, v) = 1— -. (62) 

7T [U — V) 1 + J 2 

The main assumption of the dynamical narrowing theory is that due to the small pre-factor 
in G(u>) the Lorenzian width 7 is small, so that G(u) and Q(u) are smooth on the scale 
7. With this assumption one can replace A C}S (uj, v) by 5(u — v) for v ^> 7 (because of the 
factor uj~ 2 ). Then for v > one finds 

,. M = I_^ + G(!^, ..(„) = (63) 

7T I/ 2 + 7^ ZI/ 2 2l/ 2 

On the other hand, at v ^> 7 the Lorenzian and the term with G(0) cancel each other and 
we find that in the far wings of the line the form-factor follows the spectral densities entering 
the dephasing factor, i.e. 



*.(") = I^p *"M = -i£r' ^>0,^»7. (64) 



It is important to note, that this is in fact a "perturbation theory" result and can be obtained 
expanding the dephasing factor in W(t) and H(t). Within the Lorenzian line, when v % 7, 
the estimates for the non-Lorenzian contributions in Eq.(^) can be obtained replacing in 
the denominators v 2 by 7 2 . They are negligible for small enough u, when 

G{v)-G(0),Q(v)<G{0). (65) 

If one can define a scale Au, which is the spectral width of G(u) and Q(oo), the results can 
be summarized in a simple way. The dynamical narrowing theory is valid when 7 Au. For 
detuning much smaller than the spectral width, v -C Ao>, the line form-factor is Lorenzian 
according to Eqs.(^) and (0). When the detuning is of the order or larger than the spectral 
width, v > Au, the form-factor follows the spectral density according to the perturbation 
theory result given by Eqs. (|B1D . In this domain the line is asymmetric, the red wing [y < 0) 
being enhanced compared to the blue one [y > 0). This is because of ZPFs, which enhance 
the probability for the oscillator to loss energy compared to the probability to gain it, due to 
the interaction with the thermal bath. This asymmetry is related to the perturbation theory 
result Eq . (|35|) , where the non-symmetrized current fluctuation spectral density enters. 
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VII. THE WIDTH OF JOSEPHSON LINES 



Using the results of the previous section we can find the width and the shape of the 
Josephson lines. For the line kuj one has to put 

G(u) = k 2 g(u), Q(u) = k 2 q(u), (66) 

and as a result the Lorenzian width of this line is according to Eqs.(|52|) 

7 fe = k 2 T, hT = — 
9° 

As was mentioned already, our basic assumption about the weakness of the current 
fluctuations means that the line width is small compared to the line separation, which 
means that the result given by Eqs. (|67|) is correct when 

Ik < ujj. (68) 

We will show now that this restriction, together with p ~ 1, ensures that in the vicinity of 
the Josephson frequency, i.e. at \u — kujj\ <C ujj, the dynamical narrowing theory is valid 
and the form-factor of the Josephson line is a symmetric Lorenzian with the width given 
by Eqs.flBTp. The non-Lorenzian contributions and the asymmetry of the line appears only 
at I a? — kujj\ ~ cu j, but in this domain one has to employ perturbation theory including 
non-diagonal terms with k' 7^ k. Two things have to be checked. First, that the Lorenzian 
width 7^ is small compared to the scales of the spectral densities G(w) and Q(co), given by 
Eqs.([66|). Second, the conditions Eqs. (|65|) , i.e. that the non-Lorenzian contributions to the 
form-factor are small. 

As can be seen from the results of dynamical narrowing theory, the form-factor at given 
detuning v do not depend on the behavior of the spectral densities G(u) and Q(u) at 07 > v. 
Since we are interested only in |o> — kuoj\ ujj, the screening factor n(o>) in g{uS) and q(u) 
can be neglected. As a result, it follows from Eq.(f49D that the scale of G{u) is the larger of 
T/h and ujj. The situation is more complicated for Q(u) which is linear in 07 and has no 



p 2 T H — huj coth 



2T 



(67) 
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scale. However, Eq.(|63|) for the asymmetric contribution to the form-factor is still correct 
for v 3> 7, since the effective scale of a linear Q(u) is u. 

Consider first the situation of high temperatures, T 3> hutj, when the scale of G{ui) is 
Au = T/h, and since T/h ^> uj, it follows from Eq. fl68|) , that this scale is large compared 
to 7^. Using Eq.(pE§D and Eq.(fJ^) one can check that the conditions Eq.(^) are satisfied 
when v <C T/h, which is valid automatically since v <C ujj <C T/h. In the case of low 
temperatures, T <C huij, the scale of G(uS) is Au = uj, which is large compared to 7^ 
because of Eq.fl68|). The conditions Eq. (|65|) require v <C ujj, which is also satisfied. 

The width at high temperatures, Y = (gcr 2 )^ 1 2(2p 2 + I) (T/H), was obtained in || 
considering the singular term Ak/ {oj — kuoj) 2 in (SV 2 )^ calculated from perturbation theory, 
Eq. (p5]) , as a wing of a Lorenzian line and assuming that the total energy radiated in a this 
line is not influenced by the fluctuations. In the high temperature case the fluctuations can 
be considered to be weak, Eq. fl68|) , if ga 3 ^> k 2 (T /huo). 

At low temperatures T = ((yf<T 2 ) _1 u;j, and Eq. (|68|) reduces to go 2 ^> k 2 . Finite dephasing 
of the Josephson lines exist at zero temperature due to the ZPFs, which are active because 
the junction is far from equilibrium. 

One can see from the above estimates that the fluctuations are weak due to the large 
conductance g of the shunting resistor. Approaching to the threshold p — 1, when a — > 0, the 
fluctuations increase. The effect of the fluctuations is stronger for high harmonics, k ^> 1. 

It is important to note that the general quantum result for the line width, Eq. (|67D , can 
not be obtained from the high temperature classical one simply replacing the temperature 
T by its quantum equivalent, (hu/2) coth(huj / 2T) . 

We note also that a general relation is valid, hT = (47re 2 /h)S(0), which is a simple relation 
between the width of main Josephson line at uj = ujj and the low frequency noise. As was 
mentioned in sec.|I|, there are two indications, pointing to the smallness of the fluctuations: 
the ratio of the line width to the line separation T/uj, and the ratio of the background 
contribution to the contribution of the lines, which can be estimated as S(0)u>j/(g I c R) 2 . 
It follows from the relation between T and S(0) that both ratios are of the same order, i.e. 

18 



that both ways to estimate the strength of fluctuations are equivalent. 
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